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Abstract 

In the framework of the A/r-factorization approach, the production of prompt photons 
in association with a heavy (charm or beauty) quarks at high energies is studied. The 
consideration is based on the 0(aa 2 s ) off-shell amplitudes of gluon-gluon fusion and quark- 
(anti) quark interaction subprocesses. The unintegrated parton densities in a proton are 
determined using the Kimber-Martin-Ryskin prescription. The analysis covers the total and 
differential cross sections and extends to specific angular correlations between the produced 
prompt photons and muons originating from the semileptonic decays of associated heavy 
quarks. Theoretical uncertainties of our evaluations are studied and comparison with the 
results of standard NLO pQCD calculations is performed. Our numerical predictions are 
compared with the recent experimental data taken by the D0 and CDF collaborations at the 
Tevatron. Finally, we extend our results to LHC energies. 
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Production of prompt photonqj in hadron-hadron collisions at high energies is a subject 
of pointed theoretical and experimental investigations up to now [IHZ]- Such processes 
have provided a direct probe of hard subprocess dynamics, since the produced photons are 
largely insensitive to the effects of final-state hadronization. Inclusive prompt photon cross 
sections are strongly sensitive to the parton content of a proton since, at leading order 
(LO) of Quantum Chromodynamics (QCD), such photons are produced mainly via quark- 
gluon Compton scattering or quark- ant iquark annihilation. An additional information about 
interaction can be extracted from the photon and associated jet production cross sections 
if flavor of the produced jets is tagged. Such measurements have been performed recently 
by the D0 [X i 12] and CDF J3JH] collaborations at the Tevatron where 7 + c-jet and 7 + 6- 
jet production cross sections as a function of photon transverse momentum p^ have been 
reported. The perturbative QCD predictions [7] calculated at next-to-leading order (NLO) 
agree reasonably well with the measured cross sections at relatively low p^ values, up to 
pTp ~ 70 GeV. However, the substantial disagreement between theory and data for both 
7 + 6-jet and 7 + c-jet production at large p\ was observed p~H3]. In the present note we will 
analyse recent D0 [HE] and CDF [3j|4] data as well as previous CDF measurements [HIE] using 
the /cT-factorization approach of QCD |8J. This approach has been successfully applied by us, 
in particular, to heavy flavour production E0 and inclusive prompt photon production [101413] 
at the HERA, Tevatron and LHC energies^. An additional motivation for our study is that 
such processes are background processes for the physics beyound the Standard Model (SM), 
for example the production of excited quarks or gauge-mediated supersymmetry breaking 
with neutralinos radiatively decaying to gravitinos [IH], and, therefore, it is necessary to have 
a realistic estimations of corresponding cross sections within QCD. 

First application of ^-factorization approach to production of photons associated with 
the charm or beauty quarks have been performed in our previous paper [16] . The considera- 
tion was based on the 0{aa 2 s ) amplitude for the production of a single photon associated with 
a quark pair in the fusion of two off-shell gluons g*g* — > ^QQ (see Fig. la). A reasonably 
good agreement between the numerical predictions and the Tevatron data [HUE] was obtained 
in the region of relatively low p\. where off-shell gluon fusion dominates. However, the quark- 
induced subprocesses become more important at moderate and large p\ and therefore should 
be taken into account. In the present note we extend our previous predictions [15] by in- 
cluding into the consideration two additional 0(aa1) subprocesses, namely q*q* — > •yQQ 
(Fig. lb) and q*(q*)Q —> iq(q)Q (Fig. lc), where Q is the charm or beauty quark. These 
subprocesses give a sizeble contribution in the region probed by measurements [IHI] . In 
other aspects we follow the approach described in [16J. Our main goal is to give a systematic 
analysis of available D0 and CDF data [THE] in the framework of /^-factorization. Specially 
we study different sources of theoretical uncertainties. Such calculations are performed for 
the first time. 

Let us start from a very short review of calculation steps. We describe first the evaluation 
of the off-shell q*q* — > •yQQ and q*{q*)Q — > iq(q)Q production amplitudes. We denote the 
4- momenta of incoming quarks, produced photon and outgoing quarks by ki, k 2 , p, p\ and 
P2, respectively; and polarization 4- vector of produced prompt photon is given by e. The 

1 Usually the photons are called prompt if they are coupled to the interacting quarks. 

2 A detailed description and discussion of the fey-factorization approach can be found in |14) . 
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corresponding partonic amplitudes can be written as follows: 
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Here e, e 9 and eg are the electron, light and heavy quark (fractional) electric charges, g is the 
strong charge, m q and tuq are the light and heavy quark masses, a and b are the eight-fold 
color indexes. 

When we calculate the matrix element squared, the summation over the produced photon 
polarizations is carried with J2 efl (p) e * u (p) — ~9^ U i an d the spin density matrix for all on- 
shell spinors in final state is taken in the usual form Yj u (p)u(p) = p + m. However, in the 
case of initial off-shell quarks on-shell quark spin density matrix has to be replaced by a 
more complicated expression^. To evaluate it we follow a simple approximation proposed 
in [10]. We extend the original diagram and consider the off-shell quark line as internal line 
in the extended diagram. The extended process looks like follows: the initial on-shell quark 
with 4-momentum p and mass m radiates a quantum (say, photon or gluon) and becomes 
an off-shell quark with 4-momentum k. So, for the extended diagram squared we can write: 
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3 We neglected virtualities of incoming heavy quarks in q*(q)Q —> jq{q)Q amplitude and treat them as 
on-shell. 
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where T is the rest of the original matrix element. The expression presented between T 1 ' 
and 7^ now plays the role of the off-shell quark spin density matrix. Using the standard 
on-shell condition ^Zu{p)u{p) = p + m and performing the Dirac algebra one obtains in the 
massless limit m — > 0: 

\ M \ 2 ~ (^ tr l T " ( k2 P - 2 (P ■ T J • ( 12 ) 

Now, using the Sudakov decomposition k = xp + kx (where kx is the off-shell quark non-zero 
transverse 4-momentum, k 2 = k\ = — ~k T ) and neglecting the second term in the parentheses 
in (12) in the small- a; limit, we easily obtain: 

|M| 2 ~ -^tv\np%] . (13) 

Essentially, we have neglected here the negative light-cone momentum fraction of the incom- 
ing quark. The properly normalized off-shell spin density matrix is given by xp, while the 
factor 2/xk 2 has to be attributed to the quark distribution function (determining its leading 
behavior). With this normalization, we successfully recover the on-shell collinear limit when 
k is collinear with p. 

Further calculations are straighforward and in other respects follow the standard QCD 
Feynman rulefl The evaluation of traces was done using the algebraic manipulation system 
Form [T7]. We do not list here the obvious expressions because of lack of space. The analytic 
expression for the \M\ 2 {g*g* —> jQQ) has been derived in our previous paper [T2] . 

According to the fc T -factorization theorem, to calculate the cross section of any process 
one should convolute the off-shell partonic cross sections with the corresponding unintegrated 
(dependent from the transverse momenta) quark and/or gluon distributions in a proton: 

a= J2 / ^ii( x i; x 2,k?T> k lr) fi(xu\4T^ 2 )fj( x 2,^- 2 2T^ 2 ) x 
xdx\dx2 dkfrpdklrp— — — — , 
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where cr*j(xi, X2, k^ T , k^y) ~ |A^f,(ki T , k| r ) | is the relevant off-shell partonic cross sec- 
tion and incoming partons have fractions x\ and X2 of initial protons longitudinal momenta, 
transverse momenta kiy and k 2 T and azimuthal angles <fii and 02- The multiparticle phase 
space <i$ = Ud 3 pi/2EiS^ (Y,P m ~ Y,P° nt ) is parametrized in terms of transverse momenta, 
rapidities and azimuthal angles of relevant particles: 



2Ei 2 rix * 2tt 
From (14) and (15) we can obtain following expressions 
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4 We neglected the contributions from the so-called fragmentation mechanisms. It is because after applying 
the isolation cut (see [IrEj) these contributions amount only to about 10% of the visible cross section. The 
isolation requirement and additional conditions which preserve our calculations from divergences have been 
specially discussed in [T2] . 
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where y is the rapidity of the produced prompt photon, p 1T , p 2 T, Vi, 2/2, V'i an d ^2 are 
the transverse momenta, rapidities and azimuthal angles of final state quarks, respectively. 
Similar expression can be easily written for q*(q*)Q — > iq{q)Q contribution. If we aver- 
age (16) and (17) over <fi\ and 2 and take the limit k 2 T — > and k 2T — > 0, then we recover 
corresponding expressions of usual collinear QCD approximation. 

Concerning the unintegrated quark and gluon densities in a proton, we apply the Kimber- 
Martin-Ryskin (KMR) approach [IB] to calculate them. The KMR approach is the formalism 
to construct the unintegrated parton distributions from the known conventional ones. In this 
case the unintegrated quark and gluon distributions are given by 
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where P a b(z) are the usual unregulated LO DGLAP splitting functions. The theta functions 
which appears in (18) and (19) imply the angular-ordering constraint A = |kr|) 
specifically to the last evolution step to regulate the soft gluon singularities. Numerically, for 
the input we have used LO parton densities xq(x, fi 2 ) and xg(x, /x 2 ) from recent MSTW'2008 
set 

Other essential parameters were taken as follows: renormalization and factorization scales 
/1 = ^pTp (where we vary the parameter £ between 1/2 and 2 about the default value £ = 1 
in order to estimate the scale uncertainties of our calculations), LO formula for the strong 
coupling constant a s (/i 2 ) with rif = 4 and Aqcd = 200 MeV, such that a s (M§) = 0.1232. 
We set the charm and beauty quark masses to m c = 1.5 GeV and m& = 4.75 GeV. The 
multidimensional integration has been performed by means of the Monte Carlo technique, 
using the routine VEGAS [20] • The full C++ code is available from the authors on request^]. 

We now are in a position to present our predictions. In Figs. 2 — 5 we confront the 
calculated 7 + 6-jet differential cross sections as a function of photon transverse momentum 
with the recent data [IH1] taken by the D0 and CDF collaborations at y/s = 1960 GeV. These 
data refer to the central kinematic region for heavy quark jet, and the D0 measurements [2] 
have been performed in two regions of kinematics, defined by y 7 ^ ct > and y 7 ^ ct < 0. For 
comparison we also plot the NLO QCD predictions [7] taken from [H-@]. We find that the full 
set of experimental data are reasonably well described by the fc^-factorization approach. One 
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can see that the shape and absolute normalization of measured cross sections are adequately 
reproduced. The overall agreement of our predictions and the data at moderate p\ is on 
the same level as it was achieved in the framework of NLO QCD approximation, and the 
difference appears at small and large pj> (see Figs. 2, 4 and 5). However, in the case of 
7 + c-jet production, the situation is a bit worse: we find a substantial disagreement between 
our predictions and D0 data at high pj. (see Fig. 6). As it was noted in [2], the disagreement 
can be reduced if, in particular, additional contributions involving intrinsic charm [21] are 
taken into account. The consideration of such contributions is out of our paper. Note, 
however, that very recent CDF data [3J for 7 + c-jet production are well described by the 
fc*r-factorization in a whole pt range (see Fig. 7). 

The relative contributions of g*g* — > jQQ and q*(q*)Q — > iq{q)Q are approximately 
equal to each other in the considered kinematical region. These subprocesses give a main 
contribution to the cross section up to p\ ~ 60 GeV. At higher photon transverse momenta 
the q*q* — > ^QQ subprocess dominates. The scale uncertainties of our calculations at high 
and low p\ are defined by uncertainties of corresponding leading production mechanisms 
and connected, in particular, with the unintegrated gluon and quark densities at low and 
high x. 

The CDF collaboration have reported data [5JE] which refer to the muons which originate 
from the semileptonic decays of associated charmed or beauty quarks. These measurements 
have been performed in the kinematic region defined by |?7 7 | < 0.9, pj. > 4 GeV and \r]^\ < 1 
at y^s = 1800 GeV. To produce muons from charmed and beauty quarks, we first convert 
them into a D or B hadrons using the Peterson fragmentation function [22J and then simu- 
late their semileptonic decay according to the standard electroweak theory. Additionally, the 
cascade decays b — > c — > fi have been taken into account. We set the fragmentation param- 
eters e c = 0.06 and = 0.006 and corresponding branching fractions to f(c — > fi) = 0.0969, 
f{b —>-//) = 0.1071 and f{b — )■ c —>■//) = 0.0802 [23]. Our predictions for 7 + // cross sections 
as a functions of transverse momentum pTp and azimuthal angle difference A0 (i.e. difference 
between the produced photon and muon transverse momenta) are shown in Figs. 8 and 9. 
We find that the /c^-factorization predictions slightly overestimate the measured transverse 
momentum distributions but agree with data within the uncertainties. 

Obtained perfect description of the A0 distribution is notable. The important role of such 
observables for understanding an interaction dynamics is well knowifl In particular, they 
give an additional insight into the effective contributions from higher-order QCD processes. 
In the naive LO QCD approximation, the distribution over A<p must be simply a delta 
function <5(A</> — 7r) since the produced particles are back-to-back in the transverse plane and 
are balanced in px due to momentum conservation. When higher-order QCD processes are 
considered, the presence of additional quarks and/or gluons in the final state allows the A(p 
distribution to be more spread. In the ^-factorization approach, taking into account the 
non-vanishing initial parton transverse momentum k<r leads to the violation of back-to-back 
kinematics even at leading order. Despite the fact that using the 2 — > 3 matrix elements in 
our consideration instead the 2 — > 2 ones makes the difference between the fc^-factorization 
predictions and the collinear approximation of QCD not well pronounced, the CDF data 
clearly favour the A>r-factorization results. 

6 See, for example, [14] . 
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Now we can try to extend our predictions to LHC energies. As a representative example, 
we define the kinematical region by the requirements |y 7 | < 2.5, 25 < < 400 GeV, 
\yi et \ < 2.2 and 18 < p^ < 200 GeV. Such region is coincide with one defined in recent 
analyses performed by the CMS and ATLAS collaborations [211 [25]. Our predictions for 
differential 7 + 6-jet cross sections as a function of photon transverse momentum pj, and 
rapidity y 1 are shown in Fig. 10. The corresponding experimental data at a/s = 7 TeV are 
still waited. 

To conclude, in the present note we apply the /cy-factorization approach to the anal- 
ysis of recent CDF and DO experimental data on the associated photon and heavy quark 
production. Using the off-shell (/c^-dependent) matrix elements of gluon-gluon fusion and 
quark- (anti) quark interaction subprocesses, we obtained a good agreement between our pre- 
dictions and the Tevatron data. As it was noted in p], our results agree better with the 
Tevatron data than the NLO QCD ones. It is important for further studies of small- 2 
physics, and, in particular, for searches of effects of new physics beyond the SM at modern 
hadron colliders. 
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a) b) c) 



Figure 1: Examples of Feynman diagrams contributing to associated photon and heavy quark 
production. Full set of diagrams can be obtained by permutations of quark and photon lines. 




Figure 2: The associated 7 + 6-jet cross section as a function of photon transverse momentum 
Pj, in the kinematical region defined by |y jet | < 1.5 and pip > 15 GeV at y/s = 1960 GeV. 
Left panels: the solid curve corresponds to the KMR predictions at the default scale /x = Ej<, 
whereas the upper and lower dashed curves correspond to scale variations described in the 
text. The dotted histogram represents the NLO pQCD predictions [Zj listed in pQ. Right 
panels: the different contributions to the 7 + 6-jet cross section. The dashed, dotted and 
dash-dotted curves correspond to the contributions from the g*g* — > •yQQ, q*q* — >■ jQQ and 
q*(q*)Q — > jq(q)Q subprocesses, respectively. The solid curve represents their sum. The 
experimental data are from D0 pQ. 
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Figure 3: The associated 7 + 6-jet cross section as a function of photon transverse momentum 
Pj, in the kinematical region defined by \yt\ < 1, |y jet | < 0.8 and p!p > 15 GeV at y/s = 
1960 GeV. Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions [7] 
are taken from [2]. The experimental data are from D0 [2J. 
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Figure 4: The associated 7 + 6-jet cross section as a function of photon transverse momentum 
•jp T in the kinematical region defined by \y'\ < 1, \y ict \ < 1.5 and # > 20 GeV at Vi = 
1960 GeV. Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions [7] 
are taken from [3]. The experimental data are from CDF [3]. 
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Figure 5: The associated 7 + 6-jet cross section as a function of photon transverse energy 
Ej, in the kinematical region defined by < 1.1, |r/J ct | < 1.5 and pi? > 20 GeV at 
y/s = 1960 GeV. Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions 
[7] are taken from The experimental data are from CDF jl]. 
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Figure 6: The associated 7 + c-jet cross section as a function of photon transverse momentum 
Pt in the kinematical region defined by \y">\ < 1, \y ict \ < 0.8 and # > 15 GeV at = 
1960 GeV. Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions [7] 
are taken from [7J. The experimental data are from D0 [2]. 




Figure 7: The associated 7 + c-jet cross section as a function of photon transverse momentum 
pTp in the kinematical region defined by \y^\ < 1, |y jct | < 1.5 and > 20 GeV at y/s = 
1960 GeV. Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions 
are taken from [3]. The experimental data are from CDF [3]. 
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p^GeV] P^GeV] 

Figure 8: The associated j+fi cross section as a function of photon transverse momentum jP T 
in the kinematical region defined by \rp\ < 0.9, \r]^\ < 1 and pj, > 4 GeV at y/s = 1800 GeV. 
Notation of all curves is the same as in Fig. 2. The NLO pQCD predictions [7] are taken 
from [5J. The experimental data are from CDF [5]. 
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Figure 9: The associated 7 + fi cross section as a function of azimuthal angle difference A</> 
in the kinematical region defined by |?7 7 | < 0.9, 17 < pj, < 40 GeV, \r]^\ < 1 and p^ > 4 GeV 
at y/s = 1800 GeV. Left panel: the solid and dashed curves correspond to the predictions 
of /c^-factorization approach and collinear QCD factorization based on the 0{aa 2 s ) matrix 
elements. Notation of all curves in the right panel is the same as in Fig. 2. The experimental 
data are from CDF [6]. 
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Figure 10: The associated 7+6-jet cross section as a function of photon transverse momentum 
Pt and rapidity y 1 in the kinematical region defined by \y 1 \ < 2.5, 25 < < 400 GeV, 
\yi et \ < 2.5 and 18 < p£* < 200 GeV at y/s = 7000 GeV. Notation of all curves is the same 
as in Fig. 2 (right panels). 
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